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Abstract. - We present a simple model of a vibrated box of sand, and discuss its dynamics 
in terms of two parameters reflecting static and dynamic disorder respectively. The fluidised, 
intermediate and frozen ('glassy') dynamical regimes are extensively probed by analysing the 
response of the packing fraction to steady, as well as cyclic, shaking, and indicators of the 
onset of glassy behaviour are analysed. In the 'glassy' regime, our model is exactly solvable, 
and allows for the qualitative description of ageing phenomena in terms of two characteristic 
lengths; predictions are also made about the influence of grain shape anisotropy on ageing 
behaviour. 



Vibrating sand results in very varied dynamics, ranging from glassy Jl], ||] to fluidised 
Recent experiments, e.g. have gone some way in validating the notion kI) that its 
essential features are captured by models which incorporate the fast relaxation of individual 
particles together with the cooperative rearrangements of clusters. In this Letter we present 
a simple model of a vibrated sand-box, which interpolates between the glassy and fluidised 
regimes, and is based on the generalisation of an earlier cellular automaton (CA) model 
of an avalanching sandpile. Our model shows both fast and slow dynamics in the appropriate 
regimes: in particular, it reduces to an exactly solvable model in the frozen ('glassy') regime, 
and provides one with a toy model for ageing in vibrated sand j|, ^) . 

We consider a rectangular lattice of height H and width W with N < HW grains located 
at its lattice points, shaken with vibration intensity T. Each 'grain' is a rectangle with sides 
1 and a < 1, respectively. Consider a grain (i,j) in row i, column j whose height at any given 
time is given by hij = nij^ + ariy+, with Uij- the number of vertical grains and n tj - + the 
number of horizontal grains below 
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Fig. 1 - A vertical grain needs to be tilted through the height Ah to reach the unstable equilibrium 
position and flop to the horizontal, while a horizontal grain needs to be tilted through an additional 
height AH to reach the vertical. 



• If lattice sites (i + 1, j — 1), (i + or (i + l,j + 1) are empty, grain moves there 
with a probability exp(— 1/r), in units such that the acceleration due to gravity, the 
mass of a grain, and the height of a lattice cell all equal unity. 

• If the lattice site (i — below the grain is empty, it will fall down. 

• If lattice sites (i — l,j ± 1) are empty, the grain at height h^j will fall to either lower 
neighbour, provided the height difference h^ — hi_i j±i > 2. 

• The grain flips from horizontal to vertical with probability exp(— m,ij(AH + Ah)/T), 
where rriij is the mass of the pile (consisting of grains of unit mass) above grain 

For a rectangular grain, AH — 1 — a is the height difference between the initial horizontal 
and the final vertical state of the grain. Similarly, the activation energy for a flip reads 
Ah = 6—1, where b = \/l + a 2 is the diagonal length of a grain. 

• The grain flips from vertical to horizontal with probability exp(— Ah/T). 

Hereafter we consider AH and Ah as phenomenological parameters, which will be seen to 
be related to the two main characteristic lengths in the system. 

In line with recent investigations of compaction [| [S], [HJ , we examine the behaviour of 
the packing fraction of our model, as a function of the vibration intensity T. Let N~ and N + 
be the numbers of vertical and horizontal grains in the box. The packing fraction <fi is: 

N+ - aN~ 

which we use as an order parameter reflective of the behaviour of the compactivity |l~l| . The 
vertical orientation of a grain thus wastes space proportional to 1 — a, relative to the horizontal 
one. 

We examine the response of the packing fraction for AH = 0.3, Ah = 0.05 to shaking at 
varying intensities in Fig. ^. Since the 'equilibrium' packing fraction tfioo (which we determined 
in separate runs) decreases with increasing intensity || , we plotted the difference <j> — (f>oo as a 
(logarithmic) function of time T in the figure, starting with the same initial packing fraction 
in each case. 

The dynamical response of the shaken sand-box includes three distinct regions, each illus- 
trated by representative curves in the figure: 

• a fluidised region (for r » 1), where we observe an initial increase (caused by a non- 
equilibrium and transient 'ordering' of grains in the boundary layer) of the packing 
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Fig. 2 - Plot of cj> — 4>oo versus lnt, for different values of T, indicated on the curves. Note that 
decreases with increasing F, and is thus distinct for each curve. 



fraction that quickly relaxes to the equilibrium values 0oc in each case. This over- 
shooting effect in Fig. increases with T, since grains ever deeper in the sand-box can 
now overcome their activation energy to relax to the horizontal. This inhomogeneous 
relaxation has been seen in earlier, off-lattice simulations of compaction [jl2|| . Analogous 
effects have also been observed in Ref. JDS. 

• an intermediate region (for r « 1), where the packing fraction remains approximately 
constant in the bulk, while the surface equilibrates via the fast dynamics of single-particle 
relaxation. The specific ip^ at which this occurs (0.917 here), is the single-particle 
relaxation threshold density observed in Ref. Jll|; non-equilibrium, non-ergodic, fast 
dynamics allows single particles locally to find their equilibrium configurations at this 
density. Analogous effects have been observed in recent experiments on colloids O), 
where the correlated dynamics of fast particles was seen to be responsible for most 
relaxational behaviour before the onset of the glass transition. 

• a frozen region (for r ^ 1), where the slow dynamics of the system results in a logarith- 
mic growth of packing fraction with time: 

=b(T)\nt + a, (2) 

where b(T) increases with T, in good agreement with experiment jT|. The slow dynamics 
has been identified [jl3| with a cascade process, where the free volume released by the 
relaxation of one or more grains allows for the ongoing relaxation of other grains in an 
extended neighbourhood. It includes the phenomenon of bridge collapse, which, for low 
vibration intensities, has been seen to be a major mechanism of compaction [||. As 
r decreases, the corresponding poo increases asymptotically towards the jamming limit 
0jam, identified with a dynamical phase transition in related work |fL3[ . 

We next investigate the analogue of 'annealed cooling', where T is increased and decreased 
cyclically, and the response of the packing fraction observed [[l| . The results obtained here are 
similar to those jl2| seen using more realistic models of shaken spheres, but the simplicity of 
the present model allows for greater transparency. 
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Fig. 3 - Hysteresis curves. Left: Ar = 0.1, t tap = 2000 time units. Right: Ar = 0.001, t tep = 10 5 
time units. Note the approach of the irreversibility point F* to the 'shoulder' rj am , as the ramp rate 
5r/t ta ,p is lowered. 

Starting with the sand in a fluidised state, as in the experiment [Q, we submit the sand-box 
to taps at a given intensity T for a time ttap and increase the intensity in steps of ST; at a 
certain point, the cycle is reversed, to go from higher to lower intensities. The entire process 
is then iterated twice. Figure || shows the resulting behaviour of the volume fraction <f> as 
a function of T, where an 'irreversible' branch and a 'reversible' branch of the compaction 
curve are seen, which meet at the 'irreversibility point' T* [jjj. The left- and right-hand side 
of Figure || correspond respectively to high and low values of the 'ramp rate' <5r/f tap H. As 
the ramp rate is lowered, we note that: 

• the width of the hysteresis loop in the so-called reversible branch decreases. The 
'reversible' branch is thus not reversible at all; more realistic simulations of shaken 
spheres [|J confirm the first-order, irreversible nature of the transition, which allows the 
density to attain values that are substantially higher than random close packing, and 
quite close to the crystalline limit Q . Precisely such a transition has also recently been 
observed experimentally in the compaction of rods [ p"5| . 

• the 'irreversibility point' T* approaches Fj am (the shaking intensity at which the jamming 
limit <^j am is approached), in agreement with results on other discrete models Jl6| . 

The simplicity of our model also permits us to explore the onset of 'glassy' behaviour 
between the regimes where fast and slow dynamics respectively predominate. We explore this 
via a configurational overlap function 

X {t Iei , &t) = jjJ2 0[A,j(*ref),A,j(*ref + At)]. (3) 

Here Bij(t) can take three distinct values depending on whether the lattice site at 
time t is (a) empty, (b) occupied by a + grain, (c) occupied by a — grain. We write [X, Y] = 
1 — Sx.y', i-e., 0[X, Y] = if X = Y. At is the time lag. Figure || shows results for different 
values of T, for T = 0.1,0.7,5. 

The left-hand panel (T = 0.1) shows the logarithmic behaviour characteristic of ageing, 
while the right-hand panel (T = 5) shows the quick equilibration virtually independent of 
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waiting times, which characterises the fluidised regime. The middle panel (r — 0.7) exemplihes 
the behaviour characteristic of the transition between the two regimes: there is an apparent 
'equilibration' into different metastable states depending on the waiting time t re f . 

The frozen regime is characterised by an absence of holes within the sand-box, and negligi- 
ble surface roughness. Here, our model reduces to an exactly solvable model of W independent 
columns of H noninteracting 'grains' cr n (t) = ±1, with a = +1 denoting a horizontal grain, 
and a = — 1 denoting a vertical grain. The orientation of the grain at depth n, measured from 
the top of the system, evolves according to a Markov dynamics with depth-dependent rates 

w(-l -> +1) = exp(-nA/i/r), , . 

w(+l^-l)=exp(-n(AH + Ah)/T), [> 

as rriij — n — H + 1 — i. 



The order parameter describing the mean orientation, which we hereafter refer to as 'orient - 
..Jiess', M(t) = (1 
of equation (Q) as 



edness', M(t) = (1/H) J2n=i M n (t), with M n (t) — (a n (t)), is related to the packing fraction 



v _ £±& -(1 -«) 
1 + a — (1 — a)<p 

At equilibrium, the orientedness profile is given by M„ iCq = tanh fn/(2£ eq )J , while the 
local equilibration time diverges exponentially with depth n, as r„ iCq w exp(n/^dyn)- These 
expressions involve two characteristic lengths of the model, the equilibrium length £ oq and the 
dynamical length £dyn, which read 

r r 

^ cq = AH' ^ dyn = Ah' (6) 

In the scaling regime where the height H and both lengths £ e q, £dyn are large, the mean 
orientedness is M cq « (2£ e jH) lncosh(iJ/(2^ eq )). For H < £ cq , M cq w /f/(4£ eq ) < 1: the 
system is very weakly ordered, even at equilibrium. For H ^> £ cq , M oq w 1 — (21n2)^ cq /iJ: 
the system is strongly ordered at equilibrium, except for its top skin layer, whose depth is of 
order £ oq . 
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Table I - Two different non- equilibrium regimes. 





Regime I 


Regime II 




Ceq < AO) < H 


A{t) < Ceq, H 




(ulnt > 1) 


(wlnt < 1) 


M(t) 


A(t)/H 


[A(2t)] 2 /(4^ cq ) 


S(t,s) 


1 - [A(t) - A( s )]/tf 


l-[A(2(t-a))]/2? 


C(t,s) 


1 - [2A(t) - A(t + s)]/H 


l-[A(2(t -*))]/!? 



As the equilibration time diverges exponentially with the depth, orientational order prop- 
agates down the system logarithmically slowly. More specifically, for a large but finite time t, 
only a top layer up to an 'ordering length' A(t) has equilibrated, with 

A(t)«£ dyn lnt. (7) 

We have M n (t) sa M n>e q for n < A(i), whereas M„(t) w for n > A(t). The most or- 
dered grains are situated at a depth comparable to A(i); they have a maximum orientedness 
A^ max (i) ~ tanh ((w/2) lnt), where 

_ £dyn _ _ I- a 

£cq ~ A/l ~ b - 1 

is the ratio of both characteristic lengths. 

The length £ oq is the length upto which disorder persists in the granular material when it 
has attained equilibrium. The length £dyn determines the length to which order has propagated 
in the granular material in the glassy regime. For grains where there are nearly equivalent 
orientations, even at equilibrium, there will be a large number of 'disordered' configurations 
in the top layer, since these will be almost equivalent to the strictly ordered one. (An extreme 
example would be a — ► 1, which approximately corresponds to the packing of spheres, which 
are known to be disordered even in the nominally 'equilibrium' state of random close packing.) 
The length £<jyn controls the rate at which order propagates as a function of time in the 
glassy regime of a compacted powder. Both lengths, in experimental terms, thus have the 
interpretation of the depth of the boundary layer in a vibrated granular system; in the first 
case, this description applies when equilibrium has been reached, while in the second case, 
this applies to the nonequilibrium evolution of a vibrated granular bed. 

In order for the model to exhibit interesting non-equilibrium or ageing effects, one must 
have A(t) <C H. The two-time quantities we investigate to explore ageing are the full two-time 
correlation function, S n (t,s) — (o- n (t)a n {s)), and the connected one, C n (t,s) = S n (t,s) — 
M n (t)M n (s), with < s (waiting time) < t (observation time). In terms of the overlap 
function of equation (||), we have S(t IC { + At, t le f) = l — 2x(t rc {, Ai). We are led to consider two 
different non-equilibrium regimes. In each case, the mean observables can be expressed, after 
some algebra, in terms of the ordering lengths only (see Table Q). In Regime I, the maximal 
ordering is very close to perfect, as 1 — M max (t) ~ t~ u <C 1. This is the conventional frozen 
regime (to which our data in Figure || correspond). The top layer of the system is strongly 
ordered, most of the grains are flat, and likely to stay that way: the ageing phenomenon 
corresponds to the slow ordering attempts of grains deeper in the bulk, quantified by the 
logarithmic growth of the ordering length A(t). Table || shows that the mean orientedness is 
nothing but the fraction A(t) /H of the system that has equilibrated. The two-time correlations 
are non-stationary, and they involve A(s), A(i), and A(i + s). 



(8) 
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In Regime II, the maximal ordering is very weak, as M max (i) w (w/2)lnt <C 1. This 
regime exists only for i.e., a — > 1 in the geometrical model. It corresponds to an even 

slower dynamics, since now any attempts at ordering are hindered additionally by a strong 
probability that a horizontal grain will flip to the vertical orientation. Table [j] shows that the 
mean orientedness involves the square of the ordering length, while the two-time correlations 
do not exhibit any non-stationary features characteristic of ageing, at least to leading order, 
in this scaling regime. 

The physical difference between the two scenarios is comprehensible in terms of disorder in 
grain shapes. Where grains are very irregularly shaped (Regime I), the non-equilibrium regime 
will carry all the usual characteristics of ageing. Where, however, grains are regularly shaped 
(Regime II) , the signatures of ageing will be hard to detect even in a highly non-equilibrium 
regime. It would be interesting to test this experimentally: would ageing experiments carried 
out separately on weakly vibrated rods (Regime I) or spheres (Regime II) have different 
results? 

In conclusion, the simplicity of our model makes it a useful conceptual tool for probing the 
dynamical responses of vibrated sand, from the fluidised to the frozen regimes. In the latter 
case, our model is exactly solvable, which allows one to describe the by now well-established 
picture of logarithmic compaction, in terms of two characteristic lengths. The improvement 
of this necessarily qualitative picture of ageing by the addition of more realistic and complex 
interactions, while still retaining the overall conceptual simplicity of our model, constitutes 
the focus of current research. 
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